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Cluster structure of a low-energy resonance in tetraneutron 
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We theoretically investigate the possibility for a tetraneutron to exist as a low-energy resonance 
state. We explore a microscopic model based on the assumption that the tetraneutron can be treated 
as a compound system where 3 n+n and 2 n+ 2 n coupled cluster configurations coexist. The influence 
of the Pauli principle on the kinetic energy of the relative motion of the neutron clusters is shown 
to result in their attraction. The strength of such attraction is high enough to ensure the existence 
of a low-energy resonance in the tetraneutron, provided that the oscillator length is large enough. 
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The first claim of the experimental observation of the 
nuclear stable tetraneutron has been made in Since 
then all other experimental attempts to find either a 
bound or a resonance state in the system of four neutrons 
have not met with success. However, in a recently re- 
ported experiment with the breakup of 14 Be [2j six events 
consistent with the formation of a bound tetraneutron 
were revealed. Unfortunately, several other experiments 
[1, H, Q undertaken to verify these results failed to prove 
or refute completely the existence of the tetraneutron 
due to a poor statistics. An overall conclusion of a num- 
ber of theoretical papers on this subject @, 0, H, 9 is 
that it does not seem possible to change modern nuclear 
Hamiltonians to bind a tetraneutron without destroying 
many other successful predictions of those Hamiltoni- 
ans. For instance, calculations within the hyperspherical 
functions method (HSFM) [7| suggest that a very strong 
phenomenological four-nucleon force is needed in order 
to bind the tetraneutron. And yet neither theoretical 
nor experimental results exclude the possible existence 
of tetraneutron as a low-energy resonance (see [f| HI, [l(| ) • 

There are only few cases of theoretical treatment of the 
resonant tetraneutron. In Ref. [9( the continuum states 
of the 4 n system were studied in the framework of the 
approach which combines concepts of the HSFM and the 
resonating-group method (RGM). Along with the lowest 
order hyperharmonic the authors of Ref. 0] invoked the 
hypcrharmonics with hypermomenta K = i^min + 2 and 
K = K m i n + 4, which reproduce 2 n+ 2 n clustering of the 
tetraneutron. The analysis of the energy behaviour of 
the eigenphases led authors to the conclusion that 4 n has 
a resonance state at an energy of about 1-3 MeV. But 
clear indication of such a resonance was obtained only 
for the Volkov effective NN potential, which is known to 
be inappropriate for studying multineutron systems as it 
binds dineutron. 

The most systematic study of four-neutron resonances 
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was performed in 0] , where configuration space Faddeev- 
Yakubovsky equations have been solved using realistic 
NN interaction to follow the resonance pole trajectories 
in the complex energy plane. It was concluded in [f| that 
tetraneutron - bound or resonant - can be created only in 
strong external fields and would disintegrate right after 
such a field is removed. The same authors remarked, 
however, that an accurate determination of the physical 
resonance position was not possible with the methods 
used in 

Hence, the existence of resonant tetraneutrons is still 
an open question. In the present paper, we study the 
possibility for the tetraneutron to be described as a com- 
pound system where 3 n+n and 2 n+ 2 n coupled cluster 
configurations coexist. 

The Pauli exclusion principle is known to significantly 
influence the interaction of composite particles. An exact 
treatment of the antisymmetrization effects related to the 
kinetic energy exclusively was shown to result in either an 
effective repulsion or attraction of the clusters [llj . Such 
an effective interaction substantially affects the dynam- 
ics of the cluster-cluster interaction and can, on occasion, 
produce resonance behaviour of the scattering phase or 
even a bound state in compound nuclear system. The 
role of the Pauli principle in the formation of both the 
discrete spectrum and multi-channel states of the two- 
cluster nuclear systems can be studied in the Algebraic 
Version of the resonating-group method (AVRGM) 
An RGM wave function is built in the form of an antisym- 
metrized product of cluster internal wave functions [ItJ 
and a wave function of their relative motion. The latter 
depends only on the Jacobi vector R of the considered bi- 
nary system and, according to the AVRGM, is expanded 
in a complete discrete basis of harmonic-oscillator (HO) 
states allowed by the Pauli principle 

*„ (B) (R) = ]T \/A^ B Vn(R). 

n 

The Pauli-allowed basis functions ip n are denned in the 
Fock-Bargmann representation [13j , where they take ex- 
ceptionally simple form. Expansion coefficients both of 
the discrete eigenstates with energy E K — — k 2 /2 < 



2 



and of the continuum eigenstates {C n (E)} with energy 
E > are found by solving a set of linear equations 



(n\H\n')C n > - EA n C n = 0. 



(1) 



The asymptotic behavior of the continuum eigenstates is 
expressed in terms of Hankel functions of the first and 
second kind, and the scattering S'-matrix elements. 

The set of quantum numbers n includes the number of 
oscillator quanta u, the indices (A, fi) of their SU(3) sym- 
metry, the additional quantum number otr\,n) if there are 
several differing (A, /i) multiplets, the orbital momentum 
L and its projection M, and one more additional quan- 
tum number if the multiplet (A, //) has several states 
with the same values of L. In this paper we restrict the 
discussion to the case of L = M — and, hence, we will 
use only first three quantum numbers. 

The influence of the Pauli exclusion principle on the 
collision of clusters through the kinetic energy is deter- 
mined by the eigenvalues A v of the antisymmetrizer: 

Aijj u = A v ip u , A v > 0, A u — ► 1, v — > oo. 

The eigenvalues are proportional to the probability of the 
system being in the state determined by the eigenfunc- 
tion ip v . A u > correspond to the Pauli- allowed states, 
while zero eigenvalues belong to the Pauli-forbidden 
states. A u < 1 generate repulsion of clusters, whereas 
an attraction appears in the states with the eigenvalues 
exceeding unity. 

Due to the exchange of nucleons belonging to different 
clusters, Eq. {1} resembles the Schrodinger equation for 
relative motion of two particles with highly nonlocal in- 
teraction between them. This complicates an analysis of 
the cluster-cluster interaction generated by the nucleon- 
nucleon forces. However, if the oscillator length ro is 
much less than the range bo of the NN potential, then 
the potential energy matrix in the HO representation is 
equivalent to the diagonal matrix which is a discrete ana- 
log of cluster-cluster potential in the coordinate space 



\{v,l\U\v',l)C v , tl = U eS (r v )C v>l . 



(2) 



For a central nucleon-nucleon potential having a Gaus- 
sian form we will have 



U eS (r v ) ~ U a exp <j 



A 1 + A 2 



AiA 2 



-ro 



V2i/ + 2/ + 3, 

where Uq is the strength of the Gaussian potential, r v de- 
fines the distance between clusters (composed of A\ and 
A 2 nucleons, respectively) and I is the angular momen- 
tum of the cluster relative motion. 

The validity of relation ([2]) is demonstrated in Fig. [T] 
by the example of effective 2 n- 2 n potential generated by 
the Serber NN interaction >4] (r = 0.15 fm, 6 =1.48 
fm, v = 2k). So, we deduce from §S§ that the effective 




FIG. 1: (Color online) Effective n- n potential generated by 
the Serber interaction. Solid circles: Y2y=i{^\U\k') ■ Empty 
circles: equivalent local potential U e g(rk) (see text for de- 
tails). 
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FIG. 2: (Color online) Effective 2 n- 2 n potential generated 
by the Serber interaction versus the number of quanta k at 
different values of oscillator length ro. Curves: 1 — ro = 0.15 
fm; 2 - r = 0.25 fm; 3 - r = 0.35 fm; 4 - r = 0.5 fm. 



cluster-cluster potential decreases exponentially as oscil- 
lator length tq increases (see Fig. [5]). It would appear nat- 
ural that ro should be large enough to give a reasonable 
fit to the neutron density distribution inside dineutron 
clusters, i.e. ro > 60 ■ But then we immediately arrive 
at the conclusion that effective intercluster interaction 
derived from the NN potential is of minor importance 
and can not produce a resonance in 4 n system. We have 
used the Gaussian-type potential for simplicity, but all 
the above is certain to remain valid for any short-range 
potential. 

As regards the contribution from the kinetic energy 



3 



operator (with its exchange part), in this case it is re- 
duced to the elimination of a single Pauli forbidden state 
by virtue of the fact that for the 2 n + 2 n configuration 
the eigenvalues of all the allowed states equal to unity. 
Thus, the eigenfunctions of the kinetic energy operator 
in the region of relatively small energies are similar to the 
eigenfunctions of a particle in the field of a hard core, and, 
in the region of large energies, they are identical to the 
wave functions of free motion of a particle. Hence, there 
is no grounds to believe that the tetraneutron can exist 
as a dineutron-dineutron cluster system. Behaviour of 
the phase shift of the 2 n + 2 n scattering, shown in Fig. O 
supports this inference. Figure [3] also corroborates our 
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FIG. 3: (Color online) Phase shift of the 2 n + 2 n scattering, 
ro — bo — 1.48 fm. Dashed curve: the phase obtained by in- 
corporating the kinetic energy exchange potential exclusively. 
Solid curve: phase shift calculated with the inclusion of the 
cluster-cluster potential generated by Ripka forces [l5[. 

conclusion that the 2 n- 2 n interaction generated by the 
NN forces has a minor effect on physical observables 

For each compound nuclear system, one can identify 
several cluster configurations corresponding to a certain 
set of clusters released as a result of a nuclear reaction. 
The probability for a cluster configuration to be realized 
in the Pauli-allowed basis function of a binary cluster 
system with the minimum number of quanta was shown 
to be proportional to the eigenvalue of isolated config- 
uration [16[. Maximum eigenvalues belong to the 3 n+n 
configuration of the 4 n system. They exceed unity and 
have a maximum at the minimum number of quanta that 
can be interpreted as an attraction of neutron clusters 
at small distance between them. Taking into account 
various cluster configurations of the compound nucleus 
causes the maximum eigenvalues to increase. Moreover, 
new branches of excitation appear with particularly large 
(larger than unity) eigenvalues of the allowed states that 
is an evidence of essential attraction of the clusters due 
to exchange effects. Therefore, we have good reason to 



believe that consideration of 2 n + 2 n and 3 n + n cou- 
pled cluster configurations of the 4 n system will allow us 
to solve the problem of existence of a bound state or a 
resonance in the tetraneutron. 

The 2 n + 2 n and the 3 n + n configuration coupled 
provide two branches of the two-fold degenerate SU(3) 
representation (2k, 0). The behaviour of the eigenvalues 
belonging to these branches is shown in Fig. 2J An effec- 
tive attraction reveals itself in (2k, 0)+-branch while in 
(2k, 0) _-branch an effective repulsion takes place. The 
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FIG. 4: (Color online) Eigenvalues A(a. m ) of the allowed states 
for the 4 n system in the coupled-channel approach. 

remarkable feature of these branches is that at large dis- 
tance between clusters (k — > oo) they contain the wave 
function of both cluster configurations on equal terms: 

*(2fe,o)± ->■ {i>(2k,o)( 3 n + n) ±-0(2fc,o)( 2 «+ 2 n)} . 

Returning to the analysis of intercluster interaction gen- 
erated by NN forces, we refer to Fig.[5j Figure [5]depicts 
the effective cluster-cluster potentials generated by the 
Ripka forces and calculated for isolated cluster configu- 
rations of the 4 n system. It is seen from Fig. [5] that the 
strength of 3 n-n interaction exceeds that of 2 n- 2 n inter- 
action by the difference A„ in intrinsic potential energy 
of the 3 n cluster and two dineutrons. For the Gaussian- 
type potential A u = -0.5Fi 3 (l + 2r 2 /6 2 1 )~ 3/2 , where Vi 3 
is the strength of interaction of a neutron pair in singlet 
spin state. 

The A u causes the SU(3)-branches (2k, 0)+ and 
(2k, 0)_ to remain coupled even at large k (see Fig. 
The difference in intrinsic potential energies A u of 3 n and 
two dineutrons, along with the difference in their kinetic 
energies A t , gives rise to the threshold of the tetraneu- 
tron decay into 3 n and a neutron. Although the strength 
of such energy barrier A — A u + A t is mainly determined 
by the discrepancy between intrinsic kinetic energies of 
the clusters [19|, the latter does not couple basis func- 
tions belonging to different SU(3) branches. At the same 
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FIG. 5: (Color online) Effective intercluster potential induced 



by the Ripka interaction. Solid and empty circles: 
2 n+ 2 n cluster configurations, correspondingly. 
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FIG. 6: (Color online) Effective cluster-cluster po- 
tential Uff generated by Ripka interaction for 
the 4 n system in the coupled-channel approach. 
Curves: 1 U++ = Z£ =1 <(2*,0)+|17|(2*',0) + ); 

2 Ur = E^=i<(2fc,0)_|t/|(2fc',0)-}; 3 



^ + -=E^= 1 ((2fc > 0) + |f/|(2fc',0)_>. 



time, nonzero A u implies additional repulsion, which pre- 
vents the formation of resonance in the tetraneutron. A u 
tends to diminish with increasing the oscillator length 
tq. Hence, effective attraction generated by the kinetic 
exchange terms is able to create a low-energy resonance 
in A n only at large values of ro- Figure [71 which rep- 
resents the phase shift of the elastic 2 n + 2 n scattering 
versus dimensionless energy, sustains this assumption. A 
resonance behaviour of the phase shift is observed only at 



the values of ro of the order of 10 fm. This resonance is 
very close to zero energy above the 4 n-^ 2 n+ 2 rt threshold, 
and its width is quite large in comparison with its energy. 
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FIG. 7: (Color online) Phase shift of the elastic 2 n + 2 n scat- 
tering in the coupled-channel approach, provided that A u = 
(1) and A„ / (produced by Ripka nucleon-nucleon forces 
of the range b = 1.48 fm) at r = 106 (2); r = 5b (3) and 
r = b(4). 



Hence, our results are consistent with the experimental 
data 0]. Indeed, recently the authors of the latter paper 
discussed the description of these data by means of an 
unbound-tetraneutron resonance in [Iol |. They came to 
recognize that a broad low-energy resonance in the tetra- 
neutron can be compatible with the events observed in 
2{. Moreover, theoretical ah initio calculations of Pieper 
6| also leaves room for the occurrence of a broad 4 n res- 
onance with an energy of around 2 MeV or less. 

In conclusion, we revealed that the tetraneutron has a 
good chance to exist as a compound system where 3 n+n 
and 2 n+ 2 n coupled cluster configurations coexist. The 
influence of the Pauli principle on the kinetic energy of 
the relative motion of the neutron clusters was shown to 
result in their attraction. The strength of such attraction 
is high enough to ensure the existence of a low-energy res- 
onance in the tetraneutron, provided that the oscillator 
length is large enough. It was also demonstrated that 
increasing the oscillator length results in a depression of 
the cluster-cluster potential. For this reason the results 
are not sensitive to the choice of the phenomenological 
two-body nucleon-nucleon potential. No three- and four- 
nucleon forces are employed in the calculations. 

We believe that our theoretical predictions can provide 
fresh insight into the problem of existence of resonance 
states in pure neutron systems and may help to illumi- 
nate the nature of such states and mechanism of their 
formation. 
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